function ret=epsilon(C,e0,n)
 if(C*e0>=1)
  error("Wrong C,e0");
  return;
 endif 
 ex=(3/2)^n;
 ret=(1/C)*((C*e0)^ex);
endfunction
